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INTRODUCTION 

Electronic  devices,  such  as  for  sensing  and  navigation,  are  becoming  more  common  within 
military  projectiles.  Modeling  and  simulation  may  be  useful  to  understand  the  electromagnetic 
effects  caused  by  a  spinning  metal  projectile  flying  though  a  magnetic  field  and  help  minimize 
unfavorable  aspects  of  these  effects  on  the  electronic  devices  contained  inside  the  round.  These 
effects  include  induced  (eddy)  currents  and  changes  in  magnetic  flux.  The  eddy  current  is  generated 
within  a  conductor  when  the  conductor  moves  relative  to  a  stationary  magnetic  field,  such  as  the 
earth’s  magnetic  field.  The  eddy  current  is  also  induced  within  the  conductor  when  an  alternating 
external  magnetic  field  is  applied  to  the  conductor.  The  eddy  current  produces  the  magnetic  field, 
which  opposes  the  change  in  magnetic  flux  that  induces  such  currents,  so  the  total  magnetic  flux  is 
reduced.  The  skin  effect  occurs  due  to  the  opposing  eddy  current  induced  by  the  changing  magnetic 
field  from  the  alternating  current.  The  skin  effect  is  the  tendency  of  high  frequency  current  and 
magnetic  flux  to  concentrate  near  the  outer  edge,  or  surface,  of  a  conductor,  instead  of  flowing 
uniformly  over  the  entire  cross-sectional  area  of  the  conductor,  and  then  to  decay  toward  the  center. 
The  higher  the  spin  rate,  the  more  pronounced  the  skin  effect. 

This  eddy  current  behavior  can  be  analyzed  by  using  Ampere’s  equation  (refs.  1  and  2)  and 
Bullard’s  equation  as  shown. 

Ampere’s  equation. 


<7f+7x(|)-^xS=^  m 

If  there  is  no  external  current  density,  Je  =  0,  Ampere’s  equation  can  be  reduced  to  Bullard’s 
equation  (ref.  3): 


j---Vx(rxfi)  =  0  (2) 

a  -  Conductivity,  Siemens/meter  (S/m) 

A  -  Magnetic  vector  potential,  Weber/meter  (Wb/m) 

B  -  Magnetic  flux  density,  tesla  (T) 
v  -  Velocity,  meter/second  (m/s)  or 

Q  -  Angular  velocity,  radian/second  (rad/s)  or  revolutions  per  minute  (rpm) 
po  -  Permeability  for  free  space,  47t  x  107  Henry/meter  (H/m) 
pr  -  Relative  permeability  of  the  medium 
p  -  Permeability  of  the  medium,  p  =  po  x  pr  (H/m) 
je  -  External  current  density,  Ampere/meter  (A/m) 

Both  equations  (Ampere  and  Bullard)  show  that  the  magnetic  flux  density  is  a  function  of 
velocity,  conductivity,  and  permeability,  B  =  B(v  ,<r,p).  A  closed-form  solution  (ref.  3)  and  a  finite 
element  analysis  [FEA,  (ref.  4)]  in  the  steady  state  condition  (Je  =  0,  dB/dt  =  0),  were  performed  to 
determine  the  effects  on  the  magnetic  field  distribution  inside  and  outside  a  rotating  solid  conducting 
cylinder  (hereafter  called  cylinder)  immersed  in  a  uniform  transverse  magnetic  field  by  varying  the 
spin  rate  (rpm)  and  electromagnetic  physical  properties  [conductivity  (a)  and  permeability  (p)]  of  the 
cylinder.  Similar  studies  were  conducted  by  various  authors  using  the  analytical  method  (refs.  5  and 
6)  and  finite  element  modeling  techniques  (refs.  7  and  8). 

Two-dimensional  (2D)  finite  element  time-dependent  transient  analyses  of  an  infinitely  long 
cylinder  were  modeled  using  the  ANSYS  Maxwell  program.  A  magnetometer,  which  is  an  electronic 
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sensor  used  to  measure  the  magnitude  and  direction  (angle)  of  a  magnetic  field,  was  placed  at  the 
center  of  the  cylinder’s  axis  of  rotation  for  this  analysis  and  was  modeled  as  a  nodal  point.  Two 
magnets  (NdFe35)  were  used  to  generate  a  magnetic  field  around  the  cylinder.  The  dynamic  (time 
dependent)  behaviors  of  the  eddy  current  interaction  between  the  conductor  and  magnetic  field 
excitation  were  studied.  The  magnetic  flux  density  vector  at  the  center  of  the  cylinder’s  axis  of 
rotation  was  evaluated.  The  factors  that  may  affect  the  eddy  current  interactions  and  magnetic  field 
responses  were  analyzed  by  varying  the  spin  rate  and  the  electromagnetic  properties  (conductivity 
and  permeability)  of  the  cylinder.  The  analyses  were  done  in  three  stages. 

•  Verify  the  accuracy  of  the  transient  analysis  by  comparing  the  results  to  the 
analytical  closed-form  solution  (ref.  3)  and  a  FEA  (ref.  4)  at  a  steady-state 
condition. 

•  Compare  the  magnitude  and  phase  shift  angles  of  the  magnetic  flux  density  at  the 
center  of  the  cylinder  for  three  different  conductive  materials  at  various  spin  rates 
(3600,  7200,  10800,  and  14400  rpm).  The  phase  shift  angle  is  the  variation 
between  the  initial  (0  deg,  t  =  0)  angle  and  the  angle  affected  by  the  spin  rate 
and/or  physical  materials  of  the  cylinder. 

•  Conduct  a  parametric  study  by  varying  the  permeability  and  the  spin  rate  of  the 
cylinder  for  a  specific  conductive  material. 


THE  FINITE  ELEMENT  MODEL 

The  cylinder  and  magnetic  poles  are  assumed  to  be  infinitely  long  with  uniform  cross 
sections.  The  external  magnetic  field  is  transverse  to  the  cylinder  and  there  is  no  induced  B 
(magnetic  flux  density)  in  the  direction  normal  to  the  plane  of  the  cross  section.  The  field  patterns  in 
the  entire  device  can  be  analyzed  by  modeling  the  field  patterns  in  its  cross  section.  Therefore,  a  2D 
finite  element  model  is  used  for  the  analysis. 

Figure  1  shows  the  boundary  conditions  and  magnetic  field  excitation.  The  balloon  (no 
fringing  at  infinity)  boundary  condition  is  used  for  the  model.  Two  magnetic  poles  are  used  to 
generate  a  constant  magnetic  field  for  the  cylinder  at  the  center  location.  The  moving  objects  are 
separated  from  stationary  objects.  The  center  region  of  the  geometry,  defined  as  a  band  (moving) 
object,  contains  the  cylinder  and  part  of  the  air  gap  and  rotates  relative  to  the  global  fixed  coordinate 
system  of  the  magnetic  poles.  The  band  object  is  decoupled  from  the  magnetic  poles,  and  the 
surrounding  air  and  the  regions  are  meshed  independently.  This  is  done  so  the  center  can  freely 
rotate  since  otherwise  the  mesh  definitions  would  be  violated  by  the  rotation  of  the  cylinder. 
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Boundary  Condition:  Balloon 


Boundary  Condition:  Balloon 


Figure  1 

Model  of  boundary  conditions  and  magnetic  field  excitations 

The  objective  of  Stage  1  of  the  analyses  is  to  verify  the  analysis  results  by  comparing  them  to 
the  closed-form  solution  (ref.  3)  and  FEA  (ref.  4)  with  the  same  of  permeabilities  and  spin  rates.  In 
Stage  2,  three  different  electromagnetic  properties  are  used  for  the  cylinder  to  study  the  effects  of 
the  magnetic  responses  of  the  cylinder  in  the  time  history.  Table  1  summarizes  the  electromagnetic 
properties  of  the  parts  in  the  analysis.  In  Stage  3  of  the  analysis,  how  the  permeability  and  spin  rate 
of  the  cylinder  affect  the  magnetic  responses  of  the  cylinder  is  studied.  The  cylinder  is  spun  in  a 
clockwise  rotation  in  a  range  of  3600  to  14400  rpm,  in  3600  rpm  increments. 

Table  1 

Electromagnetic  properties 


Part 

Material 

Relative 

Permeability 

Electrical 

Conductivity 

S/m 

Magnetic 

Coercivity 

A/m 

Magnet 

NdFe35 

1.09977 

6.25E5 

-8.9E5 

Conductor 

Structural  Steel 

1 

4.032E6 

0 

Conductor 

Aluminum 

1.000021 

3.8E7 

0 

Conductor 

Copper 

0.999991 

5.8E7 

0 

Vacuum 

Vacuum 

1 

0 

0 

ANALYSIS  RESULTS 

The  purpose  of  this  study  was  to  determine  the  effects  on  the  magnetic  field  distribution  in  a 
cylinder  immersed  in  a  uniform  transverse  magnetic  field  by  varying  the  spin  rates,  electrical 
conductivities,  and  magnetic  permeabilities  of  the  cylinder  in  the  time  domain.  The  magnitude  and 
the  phase  shift  angle  of  the  magnetic  flux  density  vector  at  the  center  of  the  cylinder  were 
investigated. 
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Stage  1  -  Verification  of  Finite  Element  Method 

In  order  to  verify  the  finite  element  model,  analyses  were  performed  using  these  parameters 
and  compared  to  the  closed-form  solutions  (ref.  3)  and  FEA  results  (ref.  4). 

•  The  relative  permeability  of  the  cylinder  -  one  and  1 000. 

•  The  conductivity  of  cylinder  -  3.774E7  S/m  (aluminum). 

•  Magnetic  Reynolds  number,  equation  3  -  five  and  25. 

Rm  -  R2/fijucr  (3) 

•  Angular  velocity  can  be  determined  from  equation  4. 

Q=  R2 /Rm\io  (4) 

Rm  -  Magnetic  Reynolds  number 
R  -  Radius  of  rotating  cylinder,  meter 
D,  -  Angular  velocity,  rad/s  or  rpm 

•  Spin  rate  (rpm)  and  angular  (rotational)  frequency  (Hz)  are  shown  in  table  2. 


Table  2 

Spin  rates  and  angular  frequency 


Relative 

Permeability 

Spin  Rate 
RPM 

Angular  Frequency 
Hz 

1 

6240 

104 

1 

31200 

520 

1000 

6.24 

0.104 

1000 

3.12 

0.520 

Figures  2  and  3  show  the  magnitude  of  the  magnetic  flux  density,  B,  at  the  center  of  the 
cylinder  versus  time.  The  magnetic  flux  density  reached  steady-state  condition  after  0.01  and  4.0  sec 
for  a  relative  permeability  of  one  and  1000,  respectively.  The  magnetic  flux  lines  on  the  outside  and 
inside  of  the  cylinder,  at  steady-state  condition,  closely  match  the  analytical  closed-form  solution  (ref. 
3)  and  FEA  (ref.  4)  for  all  the  cases,  as  shown  in  figures  4  and  5.  The  magnitude  of  the  magnetic 
flux  density  vector  at  the  center  of  cylinder  decreased  at  higher  spin  rates.  A  longer  time  was  needed 
for  the  magnetic  flux  density  at  the  center  of  the  cylinder  to  reach  the  steady-state  condition  when 
there  was  a  higher  permeability. 
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Figure  2 

Magnitude  of  magnetic  flux  density  at  center  of  cylinder  of  relative  permeability  1 .0 


Figure  3 

Magnitude  of  magnetic  flux  density  at  center  of  the  cylinder  of  relative  permeability  1000 
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Magnetic  flux  lines  for  the  cylinder  of  relative  permeability  1 .0 


Steady  State  Analysis 


Rm  =  0 

Transient  Analysis 


T  =  Os 

T  =  4.0s 

T  =  4.0s 

Rm  =  0 

Rm  =  5 

Rm  =  25 

Figure  5 

Magnetic  flux  lines  for  the  cylinder  of  relative  permeability  1000 
Stage  2  -  Parametrical  Sweep,  Electrical  Conductivity 

Three  different  conductive  materials  (structural  steel,  aluminum,  and  copper)  were  chosen  to 
study  the  effects  of  the  electromagnetic  responses  of  the  cylinder.  These  three  materials  have  the 
same  values  for  permeability  and  permittivity  except  they  have  different  electrical  conductivities  as 
shown  in  table  1 .  Transient  analyses  were  performed  to  determine  the  effect  of  conductivity  on  the 
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magnitude  and  the  phase  shift  angle  of  the  magnetic  flux  density  vector  at  the  center  of  the  cylinder 
and  the  magnetic  field  distribution  inside  the  cylinder  in  time  domain.  The  cylinder  was  spun  in  a 
clockwise  rotation  in  a  range  of  3600  to  14400  rpm,  in  3600  rpm  increments. 

The  magnetic  flux  density  ( B )  is  a  vector.  It  has  two  components  for  a  2D  system.  The  two 
components  are  in  the  x  and  y  directions  (B  -B^  +  By).  The  magnitude  {Bmag)o\  the  magnetic  flux 
density  and  phase  shift  angle  (0)  can  be  calculated  this  way  (eq.  5  and  6)  and  is  illustrated  in  figure 
6. 


B 


mag 


(5) 


6  =  atanliBy.Bx) 


(6) 


Phase  shift  angle  and  position  of  magnetic  flux  density  vector 

Because  in  this  analysis  the  cylinder  is  rotated  in  a  clockwise  direction,  the  phase  shift  angle 
will  also  rotate  in  a  clockwise  direction,  starting  from  quadrant  4.  As  is  shown  in  figure  6,  when  the 
phase  shift  angle  increases  from  quadrant  3  to  quadrant  2,  the  phase  angle  becomes  positive  since 
it  is  calculated  from  the  shortest  route  from  the  positive  x  axis. 

Figures  7  to  9  show  the  magnetic  flux  density  distributions  and  flux  lines  for  all  three 
conductive  materials  at  various  spinning  rates  after  reaching  steady-state  condition  (t  =  0.04  s), 
respectively.  The  maximum  magnetic  flux  density  responses  have  a  tendency  to  concentrate  at  the 
outer  surface  of  the  cylinder  with  the  highest  magnitude  in  the  copper  cylinder,  which  has  the  highest 
conductivity.  The  effect  increases  with  spinning  rate. 
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•  «  ©  © 
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■  Maximum 


Figure  7 

Magnetic  flux  density  B  distribution  and  flux  lines  for  structural  steel  cylinder  case 


Time  =  0.04s 
Speed  =  3600  rpm 

B  Distribution  Ranges 
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Minimum 
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B  Distribution  Ranges 
(0.192- 0.567)  Tesla 
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Figure  8 

Magnetic  flux  density  B  distribution  and  flux  lines  for  aluminum  cylinder  case 
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Speed  =  10800  rpm 


Time  =  0.04s 
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B  Distribution  Ranges 
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Figure  9 

Magnetic  flux  density  B  distribution  and  flux  lines  for  copper  cylinder  case 

Figures  1 0  through  1 5  show  the  magnitudes  and  phase  shift  angles  of  the  magnetic  flux 
density  in  the  time  history,  at  the  center  of  the  structural  steel,  aluminum,  and  copper  cylinder, 
respectively,  at  various  spin  rates.  Tables  3  through  5  summarize  the  magnitude,  phase  shift  angle, 
and  position  in  terms  of  coordinate  quadrant  of  the  magnetic  flux  density  vector  at  the  center  of  the 
cylinder  at  various  spin  rates  for  all  three  conductive  cylinder  materials  at  steady-state  condition  (t  = 
0.04  s).  A  comparison  of  the  magnetic  field  responses  on  the  structural  steel,  aluminum,  and  copper 
cylinder  at  various  spin  rates  are  shown  in  table  6.  Figures  16  and  17  show  the  comparison  of  the 
magnitudes  and  phase  shift  angles  of  the  magnetic  flux  density  vector  at  the  center  of  the  structural 
steel,  aluminum,  and  copper  cylinder  at  a  spin  rate  of  14400  rpm.  Of  the  three  materials  analyzed, 
the  magnitude  of  the  magnetic  flux  density  vector  at  the  center  of  the  copper  cylinder  decreased  the 
most  while  exhibiting  the  largest  phase  shift  angles.  This  is  because  out  of  the  three  materials, 
copper  has  the  highest  conductivity.  A  longer  time  was  needed  for  the  magnetic  flux  density  at  the 
center  of  the  cylinder  to  reach  the  steady-state  condition  with  a  higher  spinning  rate  and  conductivity 
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Figure  10 

Magnitude  of  magnetic  flux  density  at  center  of  structural  steel  cylinder 


: 

£  -10.00 


3600  rpm 

^  7200  rpm 

^  10800  rpm 

^  14400  rpm 

Figure  1 1 

Phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  structural  steel  cylinder 
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Magnitude  of  magnetic  flux  density  at  center  of  aluminum  cylinder 


Figure  13 

Phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  aluminum  cylinder 
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Figure  14 

Magnitude  of  magnetic  flux  density  at  center  of  copper  cylinder 


Phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  copper  cylinder 

Table  3 

Phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  structural  steel  cylinder 


Structural  Steel 

Center  of  Cylinder 

At  t  =  0.04  s 

Spin  Rate  Bx  By  Bmag  Phase_Shift  Coordinate 


(rpm) 

(Tesla) 

(Tesla) 

(Tesla) 

(deg) 

Quadra  nt 

3600 

0.380 

-0.029 

0.381 

-4.43 

4 

7200 

0.373 

-0.058 

0.377 

-8.80 

4 

10800 

0.361 

-0.084 

0.371 

-13.10 

4 

14400 

0.347 

-0.107 

0.363 

-17.22 

4 
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Table  4 

Phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  aluminum  cylinder 


Aluminum 

Center  of  Cylinder 

At  t  =  0.04  s 

Spin  Rate 
(rpm) 

Bx 

(Tesla) 

Bx 

(Tesla) 

D 

Dmag 

(Tesla) 

Phase_Shift 

(deg) 

Coordinate 

Quadrant 

3600 

0.258 

-0.212 

0.334 

-39.35 

4 

7200 

0.088 

-0.236 

0.252 

-69.53 

4 

10800 

0.005 

-0.186 

0.186 

-91.54 

4 

14400 

-0.045 

-0.131 

0.139 

-108.84 

3 

Table  5 

Phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  copper  cylinder 


Copper  Center  of  Cylinder 

At  t  =  0.045 


Spin  Rate 
(rpm) 

(Tesla) 

Bx 

(Tesla) 

Bmag- 

(Tesla) 

Phase_Shift 

(deg) 

Coordinate 

Quadrant 

3600 

0.161 

-0.245 

0.293 

-56.71 

4 

7200 

-0.013 

-0.188 

0.188 

-93.88 

3 

10800 

-0.051 

-0.102 

0.124 

-119.79 

3 

14400 

-0.066 

-0.054 

0.085 

-140.74 

3 

Table  6 

Comparison  of  magnetic  field  responses  on  the  cylinders  at  various  spinning  rates 


Structural  Steel 

Aluminum 

Copper  | 

Cylinder 

Center  of 

Cylinder 

Center  of 

Cylinder 

Center  of 

Cylinder 

Cylinder 

Cylinder 

Spin  Rate 
(rpm) 

Bmag 

(Tesla) 

Bmag 

(Tesla) 

Phase  Shift 
(deg) 

Bmag 

(Tesla) 

Bmag 

(Tesla) 

Phase  Shift 
(deg) 

Bmag 

(Tesla) 

Bmag 

(Tesla) 

Phase  Shift 
(deg) 

3600 

0.367-0.399 

0.380 

-4.43 

0.262-0.502 

0.334 

-39.35 

0.219  -  0.543 

0.293 

-56.71 

7200 

0.354-0.412 

0.377 

-8.80 

0.192-0.567 

0.252 

-69.53 

0.143  -0.602 

0.188 

-93.88 

10800 

0.341-0.425 

0.371 

-13.10 

0.145-0.599 

0.186 

-91.54 

0.096-0.625 

0.124 

-119.79 

14400 

0.328-0.438 

0.363 

-17.22 

0.111-0.615 

0.139 

-108.84 

0.068-0.637 

0.085 

-140.74 
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Figure  16 

Comparison  of  magnitude  of  magnetic  flux  density  vector  at  center  of  cylinder 


Figure  17 

Comparison  of  phase  shift  angle  of  magnetic  flux  density  vector  at  center  of  cylinder 
The  analyses  results  are  summarized: 

•  An  increase  in  conductivity  draws  the  magnetic  field  toward  the  outer  surface  of 
the  cylinder  due  to  the  skin  effect. 

•  An  increase  in  conductivity  decreases  the  magnitude  of  the  magnetic  flux  density 
vector  at  the  center  of  the  cylinder  due  to  the  skin  effect. 

•  An  increase  in  conductivity  increases  the  phase  shift  angle  of  the  magnetic  flux 
density  vector  at  the  center  of  the  cylinder  (in  this  analysis,  increasing  in  a 
clockwise  rotation  from  quadrant  4). 

•  The  effects  on  the  magnetic  field  responses  inside  the  cylinder  due  to  eddy 
currents  and  skin  effect  are  more  pronounced  when  there  is  a  higher  spin  rate. 
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•  A  longer  time  was  needed  for  the  magnetic  flux  density  at  the  center  of  the 
cylinder  to  reach  the  steady-state  condition  for  higher  spinning  rate  and 
conductivity. 

Stage  3  -  Parametric  Sweep,  Magnetic  Permeability 

The  aluminum  cylinder  was  chosen  for  this  study.  A  parametric  study  was  conducted  by 
varying  the  relative  permeability  (from  one  to  six,  in  increments  of  one)  of  the  cylinder  and  spinning  it 
in  a  clockwise  rotation  from  3600  to  14400  rpm,  in  3600  rpm  increments.  Figures  18  to  20  show  the 
resulting  plots  of  magnetic  flux  density  and  flux  lines  with  varying  relative  permeabilities  and  spin 
rates.  Figures  21  to  28  show  the  magnitudes  and  phase  shift  angles  of  the  magnetic  flux  density 
vector,  B,  at  the  center  of  the  cylinder  at  various  relative  permeabilities  and  spin  rates  in  the  time 
domain.  Tables  7  to  9  summarize  the  magnitude,  phase  shift  angle,  and  position  in  terms  of 
coordinate  quadrant  of  the  magnetic  flux  density  vector  at  the  center  of  the  cylinder  at  various 
relative  permeabilities  and  spin  rates. 
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Figure  18 

Magnetic  flux  density  B  distribution  and  flux  lines  for  relative  permeability  =  1  and  2 
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Figure  19 

Magnetic  flux  density  B  distribution  and  flux  lines  for  relative  permeability  =  3  and  4 
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Figure  20 

Magnetic  flux  density  B  distribution  and  flux  lines  for  relative  permeability  =  5  and  6 
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Figure  21 

Magnitude  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  3600  rpm) 


Phase  shift  angle  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  3600  rpm) 


Figure  23 

Magnitude  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  7200  rpm) 
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Phase  shift  angle  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  7200  rpm) 


Figure  25 

Magnitude  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  10800  rpm) 


Figure  26 

Phase  shift  angle  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  10800  rpm) 
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Figure  27 

Magnitude  of  the  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  14400  rpm) 


Figure  28 

Phase  shift  angle  of  magnetic  flux  density  at  center  of  the  cylinder  (spinning  at  14400  rpm) 

Table  7 

Magnitude  and  phase  shift  angle  at  the  surface  and  the  center  of  the  cylinder  for  relative  permeability 

1  and  2 


Relative  Permeability  =  1  Relative  Permeability  =  2 


Cylinder 

Center  of  Cylinder 

Cylinder 

Center  of  Cylinder 

Spin  Rate 

Bmag 

Bmag 

Phase  Shift 

Coordinate 

Bmag 

Bmag 

Phase  Shift 

Coordinate 

(rpm) 

(Tesla) 

(Tesla) 

(deg) 

Quadrant 

(Tesla) 

(Tesla) 

(deg) 

Quadrant 

3600 

0.262-0.502 

0.334 

-39.4 

4 

0.283-0.859 

0.386 

-62.5 

4 

7200 

0.192-0.567 

0.252 

-69.7 

4 

0.181-1.016 

0.240 

-104.5 

3 

10800 

0.145-0.599 

0.186 

-92.2 

3 

0.117-1.085 

0.153 

-135.3 

3 

14400 

0.111-0.615 

0.139 

-109.4 

3 

0.078-1.126 

0.101 

-160.3 

3 
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Table  8 

Magnitude  and  phase  shift  angle  at  the  surface  and  the  center  of  the  cylinder  for  relative  permeability 

3  and  4 


Relative  Permeability  =  3  Relative  Permeability  =  4 


Cylinder 

Center  of  Cylinder 

Cylinder 

Center  of  Cylinder 

Spin  Rate 

Bmag 

Bmag 

Phase  Shift 

Coordinate 

Bmag 

Bmag 

Phase  Shift 

Coordinate 

(rpm) 

(Tesla) 

(Tesla) 

(deg) 

Quadrant 

(Tesla) 

(Tesla) 

(deg) 

Quadrant 

3600 

0.278-1.181 

0.379 

-81.0 

4 

0.259-1.471 

0.354 

-97.2 

3 

7200 

0.152  - 1.410 

0.204 

-131.9 

3 

0.127-1.765 

0.169 

-154.7 

3 

10800 

0.088  -  1.518 

0.116 

-168.8 

3 

0.065-  1.917 

0.087 

161.9 

2 

14400 

0.054  - 1.588 

0.069 

160.2 

2 

0.036-  2.016 

0.048 

127.0 

2 

Table  9 

Magnitude  and  phase  shift  angle  at  the  surface  and  the  center  of  the  cylinder  for  relative  permeability 

5  and  6 


Relative  Permeability  =  5  Relative  Permeability  =  6 


Cylinder 

Center  of  Cylinder 

Cylinder 

Center  of  Cylinder 

Spin  Rate 
(rpm) 

Bmag 

(Tesla) 

Bmag 

(Tesla) 

Phase  Shift 
(deg) 

Coordinate 

Quadrant 

Bmag 

(Tesla) 

Bmag 

(Tesla) 

Phase  Shift 
(deg) 

Coordinate 

Quadrant 

3600 

0.236  - 1.734 

0.324 

-111.8 

3 

0.211-1.974 

0.293 

-126.2 

3 

7200 

0.103  -  2.092 

0.139 

-177.3 

3 

0.081  -  2.398 

0.114 

163.0 

2 

10800 

0.048  -  2.289 

0.066 

135.1 

2 

0.036  -  2.640 

0.050 

110.3 

2 

14400 

0.026  -  2.420 

0.033 

97.2 

2 

0.018  -  2.804 

0.023 

69.0 

1 

The  analyses  results  are  summarized: 

•  A  longer  time  was  needed  for  the  magnetic  flux  density  at  the  center  of  the 
cylinder  to  reach  the  steady-state  condition  when  there  was  a  higher  permeability. 

•  An  increase  in  permeability  increases  the  initial  magnitude  of  the  magnetic  flux 
density  and  then  decreases  the  magnitude  of  the  magnetic  flux  density  at  steady- 
state  conditions  at  the  center  of  the  cylinder. 

•  An  increase  in  permeability  increases  the  magnetic  field  on  the  outer  surface  of 
the  cylinder  due  to  the  skin  effect. 

•  An  increase  in  permeability  decreases  the  magnitude  of  the  magnetic  flux  density 
vector  at  the  center  of  the  cylinder  due  to  the  skin  effect. 

•  An  increase  in  relative  permeability  increases  the  phase  angle  of  the  magnetic  flux 
density  vector  (in  this  analysis,  the  vector  starts  in  quadrant  4  and  rotates  in  a 
clockwise  direction.  When  the  position  of  the  vector  rotates  to  the  2nd  quadrant  or 
1st  quadrant  locations,  the  phase  angle  becomes  positive,  so  as  the  phase  angle 
increases,  the  value  of  the  angle  as  calculated  from  the  positive  x-axis,  becomes 
smaller). 

•  The  effects  on  the  magnetic  field  responses  inside  the  cylinder  due  to  eddy 
currents  and  skin  effect  are  more  pronounced  with  a  higher  spinning  rate. 
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CONCLUSIONS 

The  electromagnetic  analysis  of  a  rotating  conducting  cylinder  in  a  magnetic  field  was 
evaluated  by  comparing  the  analysis  results  to  the  closed  form  solution  and  verifying  that  the 
solutions  from  both  methods  closely  match  each  other,  as  seen  by  comparing  the  resulting  magnetic 
flux  lines  on  the  cylinder.  In  addition,  a  parametric  study  was  performed  using  a  finite  element 
method  solver  to  assess  the  factors  that  affect  the  magnetic  field  distribution  on  the  cylinder  by 
varying  the  magnetic  material  properties  (conductivity  and  permeability)  and  angular  velocity.  The 
magnitude  and  phase  angle  of  the  magnetic  flux  density  vector  at  the  center  of  the  cylinder  was 
investigated. 

Bullard’s  equation  states  that  the  magnetic  flux  density  is  proportional  to  the  conductivity, 
velocity,  and  permeability.  The  results  from  the  parametric  study  show  that  the  magnitude  and  phase 
angle  of  the  magnetic  field  vector  were  affected  by  the  angular  velocity  and  the  electromagnetic 
properties  of  the  cylinder.  The  analyses  results  also  show  that  the  magnetic  field  responses  on  the 
cylinder  exhibit  the  phenomenon  of  eddy  current  and  the  skin  effect. 

•  The  eddy  currents  produce  a  magnetic  flux  which  tends  to  oppose  the  change  in 
magnetic  flux  that  induces  such  currents,  so  the  total  magnetic  flux  is  reduced. 

•  The  higher  the  conductivity,  the  larger  the  eddy  currents  are,  and  the 
larger  the  permeability  or  the  higher  the  rotational  frequency,  the 
more  pronounced  is  the  magnetic  flux  reduction. 

•  The  skin  effect  is  the  tendency  of  high  frequency  current  and  magnetic  flux  to 
concentrate  near  the  outer  edge,  or  surface,  of  a  conductor  and  then  to  decay 
toward  the  center. 

•  The  higher  the  rotational  frequency,  the  larger  the  magnetic  flux  at 
the  outer  surface  of  the  cylinder. 

•  The  higher  the  rotational  frequency,  the  lower  the  magnetic  flux  at 
the  center  of  the  cylinder. 

•  An  increase  in  permeability,  conductivity,  and  rotational  frequency  produces  a 
larger  phase  angle  of  the  magnetic  flux  density  vector  at  the  center  of  the  cylinder. 

•  An  increase  in  permeability,  conductivity,  and  rotational  frequency  increases  the 
time  needed  for  the  magnetic  flux  density  at  the  center  of  the  cylinder  to  reach  the 
steady-state  condition.  However,  the  time  needed  to  reach  the  steady-state 
condition  is  almost  instantaneous,  and  the  magnetic  field  responses  during  that 
transient  period  may  be  disregarded. 

•  The  magnetic  field  density  and  phase  shift  angle  at  the  center  of  the  cylinder 
oscillate  at  double  the  rotational  frequency.  This  doubling  of  the  frequency  was 
seen  in  the  model  because  the  model  included  two  magnetic  poles.  In  reality, 
where  an  electrical  field  is  generated  by  the  earth,  we  are  not  certain  how  the 
rotational  frequency  will  be  affected.  At  higher  permeability  and  faster  spin  rates, 
the  amplitude  of  the  oscillating  increases  and  is  more  pronounced. 

Transient  finite  element  modeling  has  been  shown  to  be  capable  of  capturing 
electromagnetic  effects  caused  by  a  spinning  cylinder  within  a  magnetic  field,  such  as  earth’s 
magnetic  field.  Changes  in  the  magnetic  field  due  to  eddy  currents  and  skin  effects  were  seen  in 
relation  to  angular  velocity  and  material  properties.  As  electronic  devices  such  as  sensors  and 
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navigation  equipment  become  more  widespread  within  spinning  projectiles,  the  ability  to  model  and 
analyze  the  magnetic  field  generated  during  flight,  due  to  eddy  currents  and  skin  effect,  can  help 
ensure  the  electronics  within  the  round  function  accurately  and  reliably. 
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